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I. INTRODUCTION 



It is well known that in Quantum Electrodynamics (QED) the infra-red (IR) divergences 
get cancelled in suitably defined cross sections by virtue of the famous Bloch-Nordseick 
theoremfl]. According to this theorem, the divergences in virtual processes get canceled 
when the contribution of real photon emission is added. It is to be noted that this can- 
cellation takes place at the level of cross-sections and not at the level of amplitudes. The 
divergences at the amplitude level arise due to inappropriate choice of initial and final states. 
In LSZ formulation, the dynamics of incoming and outgoing particles in a scattering event is 
described by the free Hamiltonian and therefore, the initial and final states used to calculate 
the transition matrix elements are taken to be the Fock states. However, in an actual ex- 
periment, due to the finite size of the detector, the charged particle can be accompanied by 
any number of photons. The Bloch-Nordseick mechanism takes into account all states with 
any number of soft photons below experimental resolution thus leading to cancellation of 
divergences. The issue of cancellation of IR divergences at the amplitude level was addressed 
by Chungj^] who showed that the divergences in matrix elements are eliminated to all orders 
in perturbation theory if one chooses the initial and final states to be charged particles with 
a suitable superposition of an infinite number of photons. Kulish and Faddeev[3| defined the 
asymptotic states by means of an asymptotic Hamiltonian. They were the first to show that 
in QED, the asymptotic Hamiltonian does not coincide with the free Hamiltonian. Kulish 
and Fadeev(KF) constructed the asymptotic Hamiltonian V as for QED thus modifying the 
asymptotic condition to introduce a new space of asymptotic states given by 

|n;±) = fi£|ra) (1) 

where Q± is the asymptotic evolution operator and \n) is the Fock state. VL± is defined by 



0± = T exp 



o 



i / V as (t)dt 



(2) 



KF further modified the definition of S matrix and showed that it is free of IR divergences. 
In a nutshell, the method of asymptotic dynamics proposed by Kulish and Fadeev replaces 
the free Hamiltonian by an asymptotic Hamiltonian which takes into account the long range 
interaction between incoming and outgoing states and can be used to construct a set of 
coherent states as the asymptotic states. The transition matrix elements formed by using 
these states are then free of infra-red divergences. 



The KF method was used by Nelson and Butler [4|-|6|] to generate a set of asymptotic 
states in the asymptotic region of perturbative Quantum Chromodynamics(pQCD). The 
asymptotic states constructed were shown to lead to cancellation of IR divergences in certain 

n 

matrix elements in lowest order in pQCD. Greco eta/[7fl constructed a coherent state approach 
for non-abelian gauge theories and showed that matrix elements between coherent states of 
definite color are finite and factorized in the fixed angle regime. KF method was applied 
to QCD by Dahmein and Steiner 8| also who showed that the leading logarithmic behavior 
of the mass shell form factor can be derived from asymptotic quark gluon part of the QCD 
Hamiltonian. 

Relevance of coherent state formalism in light-front field theory (LFFT) was first dis- 
cussed by Harindranath and Varyj^] who showed that a coherent state may be a valid 
vacuum in LFFT. Later, it was shown [lO] in the context of LF Schwinger Model that the 
physical vacuum is gauge invariant superposition of coherent states of dynamical gauge field 
zero mode. 

A coherent state formalism for LFQED was developed by one of us in Ref. H[, henceforth 
referred to possible method to deal with the true IR divergences of LFFT. These 

true IR divergences are the bona-fide divergences of equal-time field theory and appear when 
both k + and approach zero. In addition to these, there are additional IR divergences 
in LFFT called the spurious IR divergences as they are just a manifestation of ultra-violet 
divergences of equal-time theory. It was shown in X that the true IR divergences in one loop 
vertex correction are eliminated when the transition matrix element is calculated between 
the coherent states in place of Fock states. Subsequently, it was proposed[3] to use the 
coherent state basis constructed in X for calculation of Hamiltonian matrix elements in 
Discrete Light Cone Quantization(DLCQ) method of bound state calculation as a possible 
way to avoid the vanishing energy denominators and the resulting true IR divergences. The 
method was applied to obtain the light-cone Schrodinger equation for positronium using 
coherent state basis and to demonstrate the absence of Coulomb singularity therein. The 
method of asymptotic dynamics has also been applied to LFQCD 13J to obtain a set of 
coherent states and it has been shown to lead to cancellation of IR divergences appearing 
due to vanishing energy denominators (which are actually the true IR divergences) in qqg 
vertex correction at one loop level. 

KF method leads to cancellation of IR divergences in QED to all orders. However, it is 
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well known that the Bloch-Nordseick theorem does not hold in QCD and therefore, in this 
case one does not expect to construct an all order proof of cancellation of IR divergences 
along the lines of KF method. Basically, the non-cancellation of IR divergences in QCD 
stems from the fact that asymptotic states here are bound states of quarks and anti-quarks 
and therefore the asymptotic Hamiltonian to be used in KF method should contain the 



confining potential and is not just the asymptotic Hamiltonian o 
method of asymptotic dynamics was introduced by McMullan etal 



QCD. An "improved" 



14Hl6j] which takes into 



181 ] - specifically due to the absence 



account the separation of particles also. The improved method has also been discussed in 
the context of LFQED and LFQCdQ. 

In this work, we calculate fermion self energy in LFQED up to 0(e 4 ). We extend the 
analysis of X to include the instantaneous interaction also in the construction of asymptotic 
Hamiltonian and the resulting coherent state basis. We show that the true IR divergences 
in electron mass renormalization are cancelled up to 0(e 4 ) if one uses coherent state basis 
for evaluating the transition matrix elements. 

Conventionally, LF quantization is performed in light-front gauge A + = due to its 
many advantages when applied to non-abelian theory 
of ghost fields. In this work on QED, we have also used light-front gauge. There has been 
some discussion addressing the question of gauge independence of LFQED calculations in 

to Q - In a recent work ' gauge mdepe,ldence ot n perturbatlve calculation of 

electron's anomalous magnetic moment has been verified [20J. However, in this work we 
have not addressed the issue of gauge dependence of our results as we plan to do this in a 
future work. 

The cancellation of IR divergences in covariant QED has been established to all orders 
both in the amplitude as well as the cross-section approach. Various authors 2lM24| have 
addressed the issue of equivalence of covariant and Light-Front formalism of field theory 
and therefore it may be considered unnecessary to address the cancellation of IR divergence 
in LFFT. However, the divergence structure of LFFT's is different from that of covariant 
formalism and there are issues present that still need to be addressed 25j. In particular, 
it is important to differentiate between true and spurious IR divergences. As discussed in 
Ref. 11 J13I \y\ a coherent state approach in LFFT is interesting due to following reason. 
LFFT's being based on a Hamiltonian approach, a coherent state method in LFFT may 
be useful from the point of view of extracting information about the artificial confining 
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potential which is needed in LF bound state calculations 26[ . It is well known that IR 
divergences do not cancel in QCD and the reason within the coherent state formalism is 
that the asymptotic states are not the asymptotic states of QCD but are bound states of 
quarks and antiquarks. In other words, if we use appropriate Hamiltonian of bound states 
as the asymptotic Hamiltonian and develop a coherent state approach based on it, then this 
approach would lead to cancellation of IR divergences in QCD as well. It will be worthwhile 
to understand this connection between cancellation/non-cancellation of IR divergences and 
the form of asymptotic Hamiltonian. The hope is that by understanding the structure of IR 
divergences we may be able to get some insight into the form of artificial confining potential 
mentioned in Ref. 26] which can then be used to perform the bound state calculations. 



The plan of the paper is as follows: In Section II, we present the Hamiltonian of LFQED 
and calculate the 0(e 2 ) electron mass renormalization using light-cone -time-ordered per- 
turbation theory(LCTOPT) in the standard Fock basis. We demonstrate the appearance of 
true IR divergences in the form of vanishing light-cone energy denominators. In Section III, 
we obtain the form of coherent states using the method of asymptotic dynamics. In Section 
IV, we calculate dm 2 in lowest order using the coherent state basis and show that the extra 
contributions due to emission and absorption of soft photons indeed cancel the IR diver- 
gences in 5m 2 . In Section V, we calculate 5m 2 up to 0(e 4 ) in Fock basis and identify the 
IR divergences in it. In Section VI we perform the same calculation in coherent state basis 
and show the cancellation of IR divergences in this basis. Section VII contains a summary 
and discussion of our results. In Appendix A, we set the notations and conventions and 
give some useful relations. Appendix B contains some useful properties of coherent states. 
Appendix C and D contain the details of the calculation of transition matrix element in 
Fock basis and coherent state basis respectively. 



II. PRELIMINARIES 



A. Light-Front QED Hamiltonian 



The light-front QED Hamiltonian in the li ght 



of independent degrees of freedom is given by [27 . 



ront gauge (A + = 0) expressed in terms 
28] 



P~ = H = Ho + V l + V 2 + V 3 , 



(3) 
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where 

H = J d 2 ^dx~{ l -^~ d^i + \{F 12 ) 2 - ^a+d-d k a k } (4) 
is the free Hamiltonian, 

Vi = e J d 2 ^ ± dx~^^ (5) 
is the standard 0(e) three point interaction, 

V 2 = -'-e 2 J d 2 x ± dx-dy-e(x- - y~)(la kl k )(x)^ \a n ^){y) (6) 

is an 0(e 2 ) non-local effective four-point interaction corresponding to an instantaneous 
fermion exchange and 

2 r 

Vi = ~\ J d 2 *Ldx-dy-& + 0(x)\x- - y-\(h + 0(y) (?) 

is an 0(e 2 ) non-local effective four-point interaction corresponding to an instantaneous pho- 
ton exchange. V2 and V3 are drawn as four-point interactions and a hash mark is drawn on 
the line representing the instantaneous particle. £(x) and a M (x) can be expanded in terms 
of creation and annihilation operators as 

72, 



+v(p, s)e i(p+x ~- p ^d\p, s, x + )}, (8) 
a ^ x) = J / 7f=T E e ;(9)[e-^"- q ^ ) a(g,A,x + ) + e^~-^)at(g,A,x + )], 

A— 1,2 

(9) 

where 

{b(p, s),b\p', s')} = 5(p + - p' + )5 2 ( P± - p' ± )6 ss , = {d(p, s), d\p\ s')}, (10) 

[a(q,X),at( q >,\>)} = 8(q + - q>+)8 2 (q ± - q! ± )8 xx , (11) 

These relations hold at equal light-front time x + . In terms of these momentum-space oper- 
ators, the free Hamiltonian has the form 



Ho = / d 2 p ± dp" 



^ «=±i ^ A=l,2 



(12) 
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Similarly, V\ has the form 

Vx = e J d 2 x ± dx~ J [dp)[dp][dk] J2 K'MPi s')tf(p, s') + e~ ip - x v(p, s')d(p, s')) 

s,s' ,X 

x^[e~ ip ' x u(p, s)b(p, s) + e ipx v(p, s)S(p, s)]e^(k)[e~ ik ' x a(k, A) + e ik ' x a\k, A)], (13) 

where 

V-i and V3 are given by following expressions: 



? >2 f 

V 2 = -— d 2 x ± dx-dy-dl[dp][dp][dk}[dk} 



s,s',X,X' 
A, 



e il ( x -v ) 



[e ip - x u(p, s)tf(p, s) + e- ipx v(p, s)d(p, s)}f(k)[e- ikx a(k, A) + e ik ' x a\k, A)] 7 + 
f \k)[e- ip - y u(p, s')b(p, a') + e^MP, s') rft (p, s')][e~*M^ A') + e^Vfc A')], (15) 



V 3 =j f d 2 ^ L dx-dy~dl[dp][dp}[dk}[dk] &llX n [^' x u(p,s')b ] {p,s') 

*J __/__/ 



s,s',or,a' 



+ e- ip - x v{p, s')d(p, s')] 1 + [e- ip x u{p, s)b{p, s) + e'^vfo s)S(p, s)}[e ik - y u{k, a')b\k, a') 
+ e- ik - y v(k, a')d(k, a')}j + [e- ik - y u(k, a)b{k, a) + e lk - y v{k, <r)cfl(k, a)] (16) 

where y = (x + , y~, xj_). 

B. Electron Mass Renormalization in Light-Front QED 

In light-front time ordered perturbation theory, the transition matrix is given by the 
perturbative expansion 

T = V + V 1 —-V+--- (17) 

p- - H 

The electron mass shift is obtained by calculating T pp which is the matrix element of the 
above series between the initial and the final electron states \p, s) and \p, a) and it is given 



by |28| 



Sm 2 =p + J2 T PP ( 18 ) 

s 

Note that only a = s contributes, as the fermion self energy is diagonal in spin. 
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We expand T pp in powers of e 2 as 



T nn = T^+T^ 



VP 



(19) 



In general, gives the 0(e 2n ) contribution to electron self energy correction. Here, the 
initial (or final) electron momentum is 



V 



P 



p 2 ± + m 2 

2p+ 



(20) 



In particular, 0(e 2 ) correction is obtained from 
T^=T^(p,p) = {p,s\V 1 



p--H { 



■Vi\p, s) + (p,s\V 2 \p, s) 




(P,s) 



(a) 



(P, a ) {p, s) 



hi 




FIG. 1: Diagrams for 0(e 2 ) self energy correction in fock basis corresponding to T\ 



Note that 



T$=TM(p,p)=T la + Ti 
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(21) 



where T la and Tu are 0(e 2 ) contributions from standard three point vertex and the 4- 
point instantaneous vertex and are represented by the diagrams in Fig. 1(a) and Fig. 1(b) 
respectively. We are interested in the true IR divergences which arise due to vanishing energy 
denominators in TOPT 11] . It is obvious that Ti& cannot have such IR divergences as there 



are no energy denominators involved and hence this term is not required in our discussion. 
Neglecting T^, 0(e 2 ) transition matrix element contributing to fermion self energy reduces 
to 

1 



TlaiPtP) = (P,S\V 1 — 



■Vi\p, s) 



(22) 



p--H 

To calculate T\ a we insert two complete sets of states so that the above equation becomes 

2 



Tla(p,P) 



E 

spins u i=l 



Id'p'^Kip, slVM, si, k[, Ai>(pi, si, k[, Ai 



1 



"1^2; S 2' ^2' ^2) 



p--H < 
(2/ 2 ,4^^lKb,3> (23) 



Substituting for V\ from Eq. ( IT3|) and using Eqs. (fT8|) and (1A4|) we obtain 

-Wi" 11 / tfrf 4fer- R --^) <24) 

where pi = p — ki. 

Calculating the trace, using Eq. ( 1A9|) for energy denominator and taking the limit A;^ — > 0, 
kix — > 0, we finally obtain 



Note that the denominator vanishes as A;^ — >■ 0, ku_ — >■ leading to true IR divergences |ll |. 



III. INFRARED DIVERGENCES AND THE COHERENT STATE BASIS 

It was shown in X, that the true IR divergences in one loop vertex correction get cancelled 
if one uses coherent state basis in LFQED. We will prove the same result for electron mass 



now obtain the form of coherent 
for equal time theory. In X, only 



renormalization in Section IV. For that purpose, we wil 
states in light-front formalism by the method used in Ref. 
three point vertex was used to obtain the asymptotic Hamiltonian and the corresponding 
coherent state basis. We extend the formalism developed in X by obtaining asymptotic limit 
of four point instantaneous interaction also. 

The light-front time dependence of the interaction Hamiltonian is given by 

H I {x+) = Vx{x + ) + V 2 {x + ) + V 3 (x + ) 



where 11] 



i 



VL(x+) = e W d^V^'Wi^) + <="'•" ,I+ ^! ,,, (^ (I, )] (26) 
h[ (u^) are three point QED interaction vertices : 

h? = E 6t & s ') b & s ^ A )^' -OtMp. s)<* , (27) 

s,s',\ 

h { 2 ] = Y. ^ s W k > A )^' S ')7^(P, s)4 , (28) 

s,s',X 

= E d & S 'WP> s ^ S ')7 M «(P, s)^ , (29) 



= S & S '^P> s ^ k > X WP> s')-fv(p, s)4 , (30) 

s,s',\ 

and z/| is the light-front energy transferred at the vertex h\ . For example, 

i/^=p +/c - p = 31 

p + + k + 

is the energy transfer at e&y vertex. The integration measure is given by 

ifi) = i I Mfl , (32) 



p + and p_|_ being fixed at each vertex by momentum conservation. 

At asymptotic limits, non-zero contributions to V\(x + ) come from regions where — > 0. 
It is easy to see that and are always non-zero and therefore, hi and h^ do not appear 
in the asymptotic Hamiltonian. Thus, the 3-point asymptotic Hamiltonian is defined by the 
following expression [ll| 

VU* + ) = e £ / dvMe A (k)[e-^ + h?\uM) + e^* + hl(uM)} (33) 

i=l,4 ^ 

where Q&(k) is a function which takes value 1 in the asymptotic region and is zero elsewhere. 

As shown in I, for V\ we can define the asymptotic region to consist of all points in the 
phase space for which 

where AE is an energy cutoff which may be chosen to be the experimental resolution. For 
simplicity, we shall choose a frame pj_ = 0. In this frame, the above condition reduces to 

^ + ^<A, (35) 
2k+ 2p+ K ' 

where A = p + AE. 

Thus, for all the points satisfying Eq. ( )35l) . and can be approximated by zero. This 
implies that in this region, the asymptotic Hamiltonian is different from the free Hamiltonian. 
For the present purpose, i.e. in order to eliminate the true IR divergences, we find it sufficient 
to choose a subregion of the above mentioned region as the asymptotic region. We define 
this subregion to be consisting of all points (k + , kj_) satisfying: 

ki < , (36) 
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fc+ < . (37) 

This choice of the asymptotic region leads to the asymptotic interaction Hamiltonian defined 
by Eq. (E3D with 



Contribution to asymptotic Hamiltonian from the four point instantaneous interaction 
can be obtained by taking \x + \ — > oo limit in V2(x + ). V2(x + ) is given by 



V 2 (x + ) = e 



2 ± J dvf\e-^ + mv?) + ^Wllag^ (39) 
where (za- ) are 4-point instantaneous fermion exchange vertices. For example, 
h? = E S ')HP, s ) a (^' A i) a ( fc 2, A 2 )ufo S ')7"7 + 7^(P, S )e^(fci)^ 2 (fc 2 ) , 

s,s',Ai ,A2 

(2) 

One can write the remaining seven terms in a similar manner. v\ is the light front energy 
transferred at the vertex hf\ For example, in Eq. HMD 

tf =p -- K + S- r = - E^Zlg±£h (40) 



and the integration measure is 

(41) 



di/(2) = 1 f [dpljdk^dh] 



(2nf 2 J y/\ 

At asymptotic limits, non-zero contributions to V2(x + ) come from regions where v\ — > 0. 
It can be shown easily that v^ 1 and vanish when = k^ and ku_ = k2j_, while the rest 
of the six v\ s are always non-zero. Thus, the asymptotic Hamiltonian for V 2 is defined by 
the following expression 

VW* + ) = e* £ / ^(». - h)le-'^h?H4 2) ) + e'-.'"-7 i f"(4 2) )] 2(p+ ^ i+) (42) 

Similarly, V3a S (^ + ) is obtained by taking the limit |a; + | — > 00 in Vs(x + ), where 

VAx*) =*Y.J nvfV'^Vf') + e'"i"^fef"(,f')] 2(±p+ 1 ± _ +)2 (43) 

Here, fo* (f} 3 ^) are 4-point instantaneous photon exchange vertices. One can easily verify 
that i^- 3 s are always non-zero for all i. Hence, V^(x + ) is zero in the asymptotic limit and 
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does not contribute to the asymptotic Hamiltonian. 

The asymptotic states can be defined in the usual manner by 

\n: coh) = Q±\n) , 

where \n) is a Fock state and £l± are the asymptotic Moller operators defined by 



(44) 



ni = Texp\-i j [V las {x + ) + V 2as {x + )]dx 



(45) 

Carrying out the standard procedure|3| of substituting k + = 0, kj_ = in all the slowly 
varying functions of k, and carrying out the x + integration, we arrive at the following 
expression for the asymptotic states: 



Q±\n: pi) =exp 



—e 



d 2 k ± r dk + 



[f(k,X:p)a\k,X)-f*(k,X:p)a(k,X)} 



+ e 2 / dp + d 2 P± / £ 



d 2 k 1± 



dkf 



E 



d 2 k 2± 



dki 



[gx{ku h, Ai, A 2 : p)a J (k 2 , X 2 )a(ki, Ai) - g 2 {fa, k 2 , Ai, A 2 : p)a(k 2 , X 2 )a)(ki, Xi)]p(p) 



where 



f(k,X-p) 



p ■ k \ p + J V m 2 

f(k,X:p) = r(k,X:p) J 



-k' 



gi(k 1 ,k 2 ,X 1 ,X 2 : p) = 
g 2 (h, k 2 , Ai, A 2 : p) 



4pH 



p ■ ki - p ■ k 2 + ki ■ k 2 
Ap+ 



S 3 (h - k 2 ) 



p ■ k\ - p ■ k 2 - ki ■ k 2 
p(p) = J2[ b l(P) b n(p)-dl(p)d n (p)} . 

n 

Applying the operator p(p) on the Fock state, we finally obtain 

dk + 



S%h - k 2 ) 



£L±\n: pi) =exp 



E 



d 2 k ± 



- W 2 J V2k^ 



|n: pi) 
(46) 

(47) 
(48) 



(49) 
(50) 



[f(k,X,p)a\k,X) - r(k,X,p)a(k,X)} 



E 



d 2 k 1L 



dkf 



E 



d 2 k 



21 



frt >2 W/ 2 J ^kf J ^ (2k)W J 



dk~ 
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\gi{kx, k 2 , Ai, A 2 : p)a ] (k 2 , X 2 )a(k 1 , A x ) - # 2 (&i, ^2, Aj, A 2 : p)a(A; 2 , A 2 )a t (fci, Xi) 

In particular, the one fermion coherent state is given by 

d 2 k, f dk + 



\n: (51) 



\P,<7' f(p)) = ex P 



E 



E 



d 2 ku_ f dkf 



[f(k,X: p)a\k,X) - f*(k,X: p)a(k,X)} 
dkt 



E 



d 2 k 2± 



[gi(k 1: k 2 , Ai, A 2 : p)a ] (k 2 , X 2 )a(k 1} X x ) - g 2 (k 1} k 2 , Ai, A 2 : p)a(fc 2) A 2 )a f (/ci, A^ 
Some useful properties of these coherent states are listed in Appendix B. 



\p,a) (52) 



IV. ELECTRON MASS RENORMALIZATION UPTO 0(e 2 ) IN COHERENT 
STATE BASIS 



\ fci 
\ 



fci / 



(P>«) 



FIG. 2: Additional diagrams in coherent state basis for 0(e 2 ) self energy correction corresponding 
to T 2 

In the coherent state basis, 0(e 2 ) self energy contribution is given by + T'W, where 
is defined in Eq. (EU) and T'W arises from 0(e 2 ) term in 



r'(p,p) ={p,s: f(p)\Vi\p,s: f(p)} 



(53) 



|p, s: f(p)) being the coherent state given by Eq. (|52|) . The contribution of 0(e) term in 
f(p) leads to additional diagrams shown in Fig. 2(a) and Fig. 2(b) and is denoted by 
T' 2a -\-T' 2h . Here a dotted line represents the soft photon in the coherent state. The diagrams 
in Fig. 2 correspond to emission (absorption) of soft photon by the incoming(outgoing) 
fermion, but since the emitted (absorbed) photon is soft, the two particle state containing it 
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is indistinguishable from a single fermion state. Substituting for V\ from Eq. ( |T3|) we obtain, 



T'(p,p) 



(2tt) 3 / 2 



\}[dpi][dk] 



\ u {Pn s 'i)t X { k i)u(pi,Si) 



2pf v / 2pfVW 
(p,s: /(p)|& t (p 1 ,s / 1 )6(pi,si)a(A;i,A)|p,s: f(p)}5 3 (p 1 -p 1 -k 1 ) 

+ (p,s: f(p)\b t (p 1 ,s' 1 )b(pi,s 1 )a<(k 1 ,X)\p,s: f(p)}5 3 (p 1 -p 1 + k 1 ] 



(54) 



Using coherent state properties Eqs. (1B1 [) and (1B3|) from Appendix B in Eq. (|54|) we obtain 



„2 /• J2l_ /• ,7k+ 



(27T) 3 ,/ 2p+ J 2k\ 
Using Eqs. (fTSjl . (IA4j) and calculating the trace, we obtain 

e 2 /■ ... r dkf (p ■ e(fci)) 2 e A (fci) 



(5m 



d 2 k 1± 



(55) 



(56) 



(2tt) 3 y ~ j kt p- h 

where the prime indicates the correction due to additional terms in coherent state basis. 
The energy denominator in Eq. ( 156]) vanishes in the limit A:^" — ^ 0, ki^ — > thus leading to 
IR divergences. However, adding Eqs. fl25|) and fl56l) . these true IR divergences get cancelled 
and the 0(e 2 ) electron mass correction is IR divergence free. 



V. ELECTRON MASS RENORMALIZATION UPTO 0(e 4 ) IN FOCK BASIS 

We will now calculate 0(e 4 ) electron mass correction in Fock basis. Transition matrix 
element for 0(e 4 ) correction to self energy is given by 



where 



T (2) = J> + J> + J, + J , + J , 



T 3 =(p,s 
T 4 =(p, s 
T 5 ={p,s 
T 6 =(p,s 
T 7 ={p,s 



z^rVi 1 TT Vi _ 1 TT Vx\p, s) 



p - H p - H p - H 
V 2 \p, s) 



1 1 
Vi— ttVi 



p - Ho p - Ho 

Vi 1 TT v 2 - 



p - Ho p -H { 

v 2 — l —-v x 



p- - Ho p- - Ho 

V 2 l —n-V 2 \p,s) 

p - Ho 



■Vi\p, s) 
Vi\p, a) 



(57) 

(58) 
(59) 
(60) 
(61) 
(62) 
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These matrix elements correspond to Figs. 3-7 and can be evaluated in the standard manner 
by inserting appropriate number of complete sets of intermediate states. We give the details 
of the calculation in Appendix C and present here the results. 




(p. ») (p. " ) 

(<■) 

FIG. 3: Diagrams for 0(e 4 ) self energy correction in fock basis corresponding to T 3 
T 3 is given by 

T 3 = T 3 (p,p) = T 3a + T 3b + T 3c (63) 

where T 3a , T 3b and T 3c correspond to Figs. 3(a)-(c) and are given by Eqs. (IClj) . (IC4j) and 
(IC6j) respectively. Using expressions for energy denominator in Eqs. ( 1A9[) and ( lAlOj) and 
using Eq. ffT8]) we obtain 



(5m% a = - / £k 1± d^ 2± f dkt dkt 



2(2^) 6 J ^ ^ J k+ k+ 
Tr[/ Al (fc 1 )(j) 1 + m) / A2 (fe 2 )(j> 2 + m) / Aa (fc 2 )(A + m) 4 Xl (h)(^ + m)] 
?,2{p-k l f{{p-k l ) + {p-k 2 )-{k l -k 2 )} 

e4 f j2i, j2i, /" dkf dk 2 



(64) 



Tr[^ A2 (A: 2 )(^3 + m) / Al (fciK^g + m) / Aa (A: 2 )( A + m) ^(fci)(^+m)] 
32(p • h)(p ■ k 2 ) [(p • A4) + (p • fc 2 ) - (Aa ■ fc 2 )] 

Tr[/ A2 (fc 2 )(^ 3 + m) / A2 (fc 2 )(j) 2 + m) ^ Al (fci)(A + ™) ^(fci)f> + m)j 

32(p- ki){p- k 2 )(p- -p' 2 ) 
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(65) 



(66) 



where p' 2 = p and p±, P2 and p 3 have been defined in Eqs. (1C2I) . (1C3|) and (1C5j) respectively. 
Note that (5m 2 )3 a , ((^m 2 ^ and (<5m 2 )3 C can have IR divergences when 

I p ■ ki -> i.e kt -> 0, k 1± ->■ 0, but p ■ k 2 ^ 0. 

II p- k 2 ->0 i.e fc^ -> 0, k 2 _L ->■ 0, but p-k x ^Q. 

III p • fci ->■ and p • fc 2 ->■ i.e. ->■ 0, ki_L ->■ 0, k 2 ->■ 0, k 2 ± -» 0. 

Now we will consider the contribution of Figs. 3(a)-(c) in each of these limits. 

Case I : In the limit kf — > 0, kxj_ — > (but p- k 2 ft 0), the contribution to T 3 from diagrams 

in Figs. 3(a) and 3(b) is given by 



[(8m 2 ) 3a + (8m 2 



1 36 1 



d 2 ki^d 2 k 2 j_ 



dA;+ dkt [2(p ■ e(k 1 )) 2 (p • e(k 2 )) 2 - (p • fc 2 )(p ■ e(/c 1 )) 2 ] 



and contribution from Fig. 3(c) is given by 



4{p ■ h) 2 {p ■ k 2 ) 



(67) 



[(8m 2 



1 3d 



(27T) 



d 2 k 1± d 2 k 2± 



dA;+ ciA; 2 + 2(p ■ e(A; 1 )) 2 (p • e(k 2 )) 2 - (p ■ k 2 )(p ■ e(k 1 )) r ' 



P 



+ 



Pt 



8^ 



.(p-hy( P -k 2 ) • (p-k^p-hyi (68) 

Here we have used Heitler method 29[ illustrated in Appendix C to deal with the vanishing 
denominator (p~ — p' 2 ~). 

Case II : In the limit k 2 — > 0, k 2 ^ — > (but p-k\ ft 0), the contribution to T3 from diagrams 
in Figs. 3(a) and 3(b) is given by 



[(8m 2 ) 3a + (8m 2 ) 3b ] H 



d 2 ki^a? 2 k 2 _i_ 



dk+ dk+ [2(p ■ e(k l )f(p ■ e(k 2 )) 2 - (p ■ h)(p ■ e(k 2 ))' 



(2tt) 6 J ~ J kt kt 

and contribution from Fig. 3(c) is given by 



4(p- h)(p- k 2 y 



(69) 



[(8m 2 " 1/7 



)3c\ 



(2n) 



d 2 ki±d 2 k 2 j_ 



dkt dkt 2(p ■ e(k 1 )f(p ■ e(k 2 )f - (p ■ h)(p ■ e(k 2 )f 



Pi 



.+ 
•1 



ip^ 



(p ■ h) 2 (p ■ k 2 ) ' (p ■ h)(p ■ k 2 ) 2 \ ^ 
Case III : In the limit kt — > 0, ku_ — > 0, kt — > 0, k 2 j_ —> 0, the contribution to T 3 from 
diagrams in Figs. 3(a) and 3(b) is given by 



[(8m 2 ) 3a + (8m 2 ) 



3/j 



ill I 



(27T) 



rf 2 ku_d 2 k 2 j_ 



dkt dkt (p-e(k 1 )) 2 (p-e(k 2 )f 
kt kt 2(p-k 1 ) 2 (p-k 2 ) 



(71) 
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and contribution from Fig. 3(c) is given by 



1 3c 



d 2 k\±d 2 k 2 ± 



(2tt) 6 y ~ j k+ k+ 

{p ■ eik^fip ■ e{k 2 )f- 



(p ■ eik^jp ■ e(k 2 )f 
4(p-A;i) 2 (p.Afc) 



+ 



4(p- fcO(p- fc 2 ) 2 



(72) 




(W 

*1 tl 



(r) 



FIG. 4: Diagrams for 0(e 4 ) self energy correction in fock basis corresponding to T4 

The contributions to T4, T5 and Tg come from diagrams in Figs. 4, 5 and 6 respectively. 
One can do similar calculation for the three cases by carefully taking the appropriate limits 
of corresponding expressions. Below we give contributions from these diagrams in each of 
the three limits. 

Case I : In the limit kf — > 0, ku_ — > (but p ■ k 2 -ft 0), the contributions of diagrams in 
Figs. 4(a), 5(a), 5(b) and 6(a) are given by 



[(Sm 2 ) ia Y = [(5m 2 ) ea y 



d 2 k 1± d 2 k 2± 1 dk ' dk2 



(2tt) 6 J J kt k+p+ 

[2p+(p ■ e(h))(p ■ e(k 2 )Mh) ■ e(k 2 )) + fc+(p ■ e(h)) 2 } 
8(p- k x )(p- k 2 ) 



(73) 



m 



I5aj 



d 2 k 1± d 2 k 2± 



dk~^ dk<2 



(2tt) 6 J " J k+ k+p^ 

[2p+(p ■ ejk^jp ■ ejk^Mh) ■ e(k 2 )) + kj{p ■ ejk^f] 
8(p- k x ){p- k 2 ) 



(74) 
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FIG. 5: Diagrams for 0(e 4 ) self energy correction in fock basis corresponding to T5 




(6) 



fa 





(p,<0 



FIG. 6: Diagrams for 0(e 4 ) self energy correction in fock basis corresponding to Tq 



[(Sm 2 ) 5b Y 



d 2 k 1± d 2 k 2± 



(2tt) 6 J " "~ J kt k+p+ 

[2pt(p ■ e(h))( P ■ e(k 2 )Mh) ■ e(k 2 )) - k+(p ■ 6(h)) 2 } 



8(p • h)(p ■ k 2 ) 

Diagrams in Figs. 4(c), 5(c) and 6(c) have IR divergences only in limit I and lead to 

dkt dk+ (p • e{h)) 2 



{(5m% c + (Sm 2 ) 6c Y 



(27T) 



d z k 1± d z k 2L 



kf k+ 2{p-k x 
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r/r 2\ U e4 /" i2i j2i f dk± dk 2 (p ■ e(kl)) , . 

i( * n w = y d k - rf k - y ^ if < 77) 

where again we have used Heitler method for evaluating T 4c + T 6c . Figs. 4(b) and 6(b) do 
not have IR divergences in this limit. 

Case II : In the limit k 2 — > 0,k 2 j_ — >■ (but p • k\ ■/> 0), the contributions of diagrams in 
Figs. 4(b), 5(a), 5(b) and 6(b) are given by following expressions 

e4 f j2,, j2i. f dkf dk£ 



i(5m 2 ), b } n = [( 5m\ b } H = - J d 2 k 1± d 2 k 2± J 



[2p+(p ■ ejk^jP ■ <k2)Mh) ■ e(k 2 )) - kj(p ■ e(fc 2 )) 5 
8(p • fci)(p • k 2 ) 

&i f J2i, j2i, /" ^1 ^2 



(78) 



[(Sm 2 ) 5a } n = ~J^J d 2 k 1L d 2 k 2L | 



fcl k 2 Pi 

[2p+(p ■ ejk^jp ■ e(fc 2 ))(e(fc 1 ) • e(fc 2 )) + fc+(p ■ e(/c 2 )) ; 
8(p • • fc 2 ) 



(79) 



[(5m 2 ) 56 ] 7J = " 7^3* / d 2 k 1± d 2 k 2± J 



k{ kjp{ 

[2pf(p ■ e(h))(p ■ e(k 2 ))(e(h) ■ e(k 2 )) - k+(p ■ e(fc 2 )) : 



(80) 



8(p ■ h)(p ■ k 2 ) 

Figs. 4(a) and 6(a) do not have IR divergences in this limit. 

Case III : In the limit kf — > 0,ku_ — > 0,k 2 — >■ 0,k 2 j_ — > 0, the sum of contributions 
corresponding to Figs. 4-6 (a) and (b) is given by 

i(5m\ a + (5m\ b + (5m\ a + (5m\ b + (5m\ a + (6m\ b ] in 

= -Wf] d k - J J4I4 (P-h)( P -k 2 ) (81) 

where we have used Eqs. flC22]) . f lC27j) . ( !C30|) and flC34l) . The last term T 7 in Eq. §2$ is 
IR convergent as the 4-point energy denominator involved here is non-zero and hence T-j is 
not needed for our discussion. One can notice that Eqs. (IS7|) - (18"T|) have true IR divergences 
since the denominator vanishes as kf — > 0, ku_ — > and/or k 2 — > 0, k 2 j_ — > 0. 



VI. ELECTRON MASS RENORMALIZATION IN COHERENT STATE BASIS 
UPTO 0(e 4 ) 

In this section, we will use coherent state basis to calculate the 0(e 4 ) electron mass 
correction. We will show that the IR divergences in additional diagrams appearing due to 
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(p. ') (p, ") 



FIG. 7: Diagrams for 0(e 4 ) self energy correction in fock basis corresponding to Tj 

the use of coherent state basis exactly cancel the IR divergences arising due to vanishing 
energy denominators calculated in Section V. In coherent state basis, 0(e 4 ) correction to 
self energy is given by 

T (2) _|_ _|_ JV _|_ T / q + J,^ 

where T 8 is 0(e 4 ) term in (p, s: /(p) | V x p _ * gp V x p _ 4 Hq V\ \p, s : f(p)) represented by Fig. 8, % 
is 0(e 4 ) term in (p, s: f (p) \ V\ p _ 4 Hq V\\p, s: f(p)) represented by Fig. 9, T{ is 0(e 4 ) term in 
(p,s: f(p)\V l —^V 2 \p J s: f(p)) + (p,s: f(p)\V 2 -=^V 1 \p, s: f(p)) represented by Fig. 10 
and T[ x is 0(e 4 ) term in (p,s: f(p)\V 2 \p,s: f(p)) represented by Fig. 11. We present the 
details of calculation in Appendix D and give below only the result for (5m 2 )'. 
The contribution corresponding to Fig. 8 is given by 

(5m% = (5m% a + (5m% b + (5m% c + (5m% d + (5m 2 )' 8e + (8m% (82) 

where {Sm 2 ) 8a - (5m 2 ) 8 j- have been evaluated in Appendix D. 
Figs. 8a and 8b contribute 

(W) -=(2^/^ 2k -/§? 

{2{p.e{k 1 )) 2 {p.e{k 2 )) 2 -\p.l 2 ){p.e{k 1 )) 2 ] 

4(p • h) 2 [(p ■ h) + (p ■ k 2 ) - (h ■ k 2 )} UA[kl) ^ 
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(a) 



fa 



(?,") (M (/) fp,o) 



FIG. 8: Additional diagrams in coherent state basis for 0(e 4 ) self energy correction corresponding 
to T 8 



86 



d 2 ki_i_d 2 k2_L 



[2(p • e(fcO) a (p • e(fc 2 )) 2 ] 



■e A (fci) 



(84) 



4(p • fci)(p • fc 2 )[(p • fci) + (p • fc 3 ) - (An • fc 2 )] 

One can notice that for (<5m 2 )§ a and (5m 2 )' 8b we need not discuss limit II as p ■ k\ is always 
small. Adding Eqs. (!83|) and ([81]) and taking the limit I, we obtain 

o4 r r jl+ ju+ 



[{5m% a + (dm 2 ) 



1 d 2 k 1± d 2 k 2± 1 dk ' dk2 



(2tt) 6 J " J kf k+ 

(p-e(h)) 2 (p.e(k 2 )) 2 (p-ejh)) 2 
2{p ■ k^ip ■ k 2 ) Aip-h) 2 



e A (fci) 



(85) 



&>,») 



fa / 



(p,«) (p.") 

(4) 

/ \ 
fc. / V 

/ \ fa 



\ 



\ 



(p,o) (P,< 



FIG. 9: Additional diagrams in coherent state basis for 0(e 4 ) self energy correction corresponding 
to T 9 
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Adding Eqs. (|67J) and §5$ we find that [(5m 2 ) 3a + (Sm 2 ) 3b Y ' +[(5m 2 )' 8a + (5m 2 )' 8b \ I is IR 
finite. 

Adding all the contributions coming from Figs. 3, 8 and 9, we find that the IR divergences 
completely cancel. Below we summarize this result for the reader's convenience: 



1. 
2. 
3. 
4. 
5. 
6. 



(5m 2 ) 3a + (5m 2 ) 3b Y+[(5m 2 y 8a + (5m 2 )' 8b Y is IR finite. 
(5m 2 ) 3a + (Sm 2 ) 3b \ n +[(5m 2 y 8c + (5m 2 )' 8d } n is IR finite. 
{5m 2 ) 3c } 1 '+[{dm 2 )' 8e } 1 is IR finite. 
{6m 2 ) 3c ] II +[{5m 2 y 8f } 11 is IR finite. 

(5m 2 ) 3 ] in + [(5m 2 )' 8a + (5m 2 y 8b ] in + [(5m 2 y 8e } in is IR finite. 
{5m 2 y 8c + (Sm 2 y 8d } in + [{5m 2 y 8f } in + [(5m 2 )' 9 ] m is IR finite. 



Thus, we can see that the self energy correction corresponding to 3-point QED vertices upto 
0(e 4 ) is IR finite. In the same manner we can show the cancellation of IR divergences for 
diagrams containing 4-point instantaneous fermion exchange vertex in all the three limits. 
We give here one calculation for illustration. The contributions corresponding to Fig. 10 is 
given by 



(8m 2 )' w = (5m 2 )' 10a + (5m 2 )' 10b + (dm 2 )' l0c + (8m 2 )' 10d + (8m 2 )' 10e 
+ (5m 2 )' 10f + (5m 2 y i0g + (Sm 2 )' 10h + (Sm 2 )' m + {5m 2 )' l0j 



As shown in Appendix D, the IR divergent contributions in Fig. 10(a), 10(c), 10(e) and 
10(g) in the limit I are given by the following expressions: 



[($rn 2 y i0a y = i(6m 2 y i0g Y =^-y 6 J d 2 k 1± d 2 k 2± J 

-Oa(^i) 
(87) 



h 'l h " 2 

[2p+(p ■ e(fci))(p ■ e(k 2 ))(e(h) ■ e(k 2 )) + fc+(p ■ e(h)) 2 ^ 
8p + (p ■ ki)(p ■ k 2 ) 



e f j2i, jai, f dk i dk z 

1 2 

[2p + (p ■ e(h))(p ■ e(k 2 ))(e(h) • e(fc 2 )) + k+(p • e(^)) 21 



[ WioJ' J d 2 k l± d 2 k 2± J 

-G>a(&i 



Bp+(p ■ ki)(p ■ k 2 ) 
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(p,s) (p,a) (p,s) (p,<r) 

(i) 



FIG. 10: Additional diagrams in coherent state basis for 0(e ) self energy correction corresponding 
to Tio 



2y 



'10eJ 



d 2 k 1± d 2 k, 



(27r)6 J - " 2± 7 A;+ 

[2pt(p ■ e(h))(P ■ <h)Mki) ■ e(k 2 )) - k+(p • e(^)) 2 ] 



dk~^ dk 2 



8pz(p ■ ki)(p ■ k 2 ) 



e A (fci) 



(89) 



Adding Eqs. (J73J175]) and Eqs. ([87] -EH]), we find that 
[(5m 2 ) 4 J'+[(5m 2 ) 5 ^+[(5m%]^^ 
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fcl / 
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(P,«) 


(P> CT ) 



fc 2 / 
/ 

/ 

/ 

/ 


\ 

\ fcl 
\ 

\ 

\ 


(P>«) 





(<■■) 



(>/) 



FIG. 11: Additional diagrams in coherent state basis for 0(e) self energy correction corresponding 
to Tn 

+ [{Sm 2 y i0g Y is IR finite. Adding all the contributions from Figs. 4, 5, 6, 10 and 11, we find 
that the IR divergences exactly cancel. Below we summarize our results: 

1. [(5m%] 7 +[(<5m%]'+[(5m^ 
+[(5m 2 )' 10g ] 1 is IR finite. 

2. [(5m 2 ) 4b ] // +[(5m 2 ) 5a ] /7 +[(5m 2 ) 5fe ]^+[(5m 2 ) 6b ]^+[(5m 2 )' 10b ]^ 
+ [(^ 2 ) / ioJ // +[(^ 2 ) / io/ 7 +[(^ 2 )W] // ^ IR finite. 

3. [(5m% c Y + i(5m% c Y + [(dm 2 )'^ + [(5m 2 )' Wl Y + [(<W%/ 



10eJ 



[(MU' + \(5m 2 )'^Y is IR finite. 
4. [(5m 2 ) 4 + (5m 2 ) 5 + (5m 2 ) 6 + (5m 2 )' 10 + (5m 2 )' n ] m is IR finite. 
Thus, finally we obtain 



(W) (2) + 5>m 2 ); 



11 



is IR finite, where (5m 2 ) ^ is 0(e A ) electron mass correction in Fock basis. This completes 
the proof of cancellation of true IR divergences upto 0(e 4 ) for fermion self energy correction 
in coherent state basis. 
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VII. CONCLUSION 



We have calculated electron self energy correction in light-front QED up to 0(e A ) and 
have shown that the true IR divergences get cancelled when coherent state basis is used to 
calculate the matrix elements. The cancellation of IR divergences between real and virtual 
processes is known to hold in equal-time QED to all orders. This cancellation was also shown 



by Kulish and Faddeev js] using the coherent state formalism. It would be interesting to 
verify this all order cancellation in LFQED. The present work is a first step towards this 
aim. The true IR divergences in QCD do not cancel in higher orders. This fact can possibly 
be used to obtain a form of the artificial potential required for the bound-state calculation. 
The connection between the asymptotic dynamics and cancellation/non-cancellation of IR 
divergences can also be exploited to explore the possibility of constructing an artificial 
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potential which is used in the bound state calculations in LFQCD 
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Appendix A: Notation and Useful Relations 
We define four-vector x M by 



The light-front variables are defined by 



x + = [ t, } , x~= 1 / \ ^ = {x\x 2 ) (Al) 



y/2 ' V2 ' 

Thus, in light-front variables 



X M = ( x +^ x j X J 
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The metric tensor is 



LLV 



10 
10 
0-10 
-1 



1. Dirac spinors: 



u(p, s) and u(p, s) satisfy the usual properties 



(ft — m)u(p, s) = 0, {t> + m)v(p, s) — 



(A2) 



u(p, s)u(p, s') = -v(p, s)v(p, s') = 2m8 ss > , u(p, s)^u(p, s') = v(p, s)^v(p, s') = 2p M <5^'(A3) 



u(p,s)u(p,s) =t> + m , v (p, s ) v (p> s ) 

s=±l/2 s=±l/2 



— m . 



(A4) 



2. Photon polarizations: 



Photon polarization tensor satisfies 



5f, + p u + 5 V +Pn 



A=l,2 



V 



Some useful properties satisfied by d a p(p) are 
tV«Mp)=-2, 



(A5) 



(A6) 
(AT) 

(A8) 
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3. Energy denominators: 



We will need the following expressions for energy denominators: 

(g • h) 
p + — kf 

p ■ ki + p ■ k 2 - h ■ k 2 
p + — kf — k 2 
p ■ k± — p ■ k 2 — k\ ■ k 2 
p+ + kf - k% 
p ■ ki - p ■ k 2 + h ■ k 2 
p + — k^ + k 2 
p ■ ki + p ■ k 2 + k\ ■ k 2 
p+ + fcT" + k+ 

Appendix B: Properties of Coherent States 







P 


- K - (P-h) 


p~ 


-K 


— k 2 


- (p - ki - k 2 ) 


p- 


+ K 


— k 2 


-{p + ki- k 2 ) 


p- 


-K 


+ k 2 


- (p - h + k 2 ) 


p- 


+ ^ 


+ K 


-(p + h + k 2 ) 



(A9) 



(A10) 



The coherent state containing a fermion and superposition of infinite number of soft 
photons is denoted by |1: p^ and is defined by Eq. (1521) . Similarly, the coherent state 
containing a fermion and a hard photon is denoted by |2 : p i} ki). The coherent states |1 : p^ 



are the eigenstates of a(k, A) [ll] 



«(M|l: P ,) = - ( J )3 /% R) |l: ft ). (Bl) 



Also, 

a(M|2: Pi,h) = - (27r e )3/2 /( ^ r ^ ) |2: Pi , h) + 5 3 (k - kJS^l : Pi ) , (B2) 



and 

a\k,p)\l: Vi ) = ^ ri ^ Pl) |1 -Pi) + 12: PiM) , (B3) 

at(fc,p)|2: Pi, k) = * f {k ^ l \ l: Pi , h) + |3: Pi , k h kj) . (B4) 
Coherent states satisfy the following orthonormalization properties: 

(1: Pf,(T f \l: Pi,(Ti) = 5 {3) (pi- Pf)5 a . a (B5) 



(1: p f ,<Tf\2: Pi ,ai,ki,Xi) = ^^^^=^<5 (3) fe - p f )S aiaf (B6) 
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Appendix C: Transition matrix element in fock basis for self energy upto 0(e 4 ) 



We now calculate T 3 which is defined by Eq. (158|) and corresponds to Fig. 3. Inserting 
complete sets of intermediate states in T 3 , we obtain 

T 3 (p, p) = T 3a + T 3b + T 3c 

where 

e 4 f d 2 ki±d 2 k 2 ± f dkfdk 



(2nyj 2p+ J 32kfktptptpt 
u(p,g)[^ Al (A;i)(A + m) jt X2 {k 2 )(£> 2 + m) ^{k 2 )( A + m) / Al (k^ujp, s) 

(p- -Pi - K)(p- -P2 - K - l%)(p~ - Pi - K) 

with 



(CI) 



Pi=p-k u (C2) 
p 2 =p-k 1 - k 2 . (C3) 



Similarly, 

e 4 f (i 2 k 1 _|_d 2 k 2 _L /" dk±dk 2 

1 r - 



36 (27r)6y 2p+ J 32kfk+ptp+pj 

u(p,s)[4 X2 (k 2 )(jf) 3 + m) 4 Xl (h)(h + ™) / A2 (fc 2 )(A + m) ^(k^ujp^s) 

(p- - ^ - fcr)(p~ - Ps - K)(p~ - P2 - K - K) 

with 



(C4) 



T- 



3c 



p 3=p -k 2 (C5) 
e 4 /" d 2 k 1± d 2 k 2± /• dA;+dfc+ 

(C6) 



(2tt) 6 y 2p+ y 32/fc+fc+p+p+p+ 

u{p,s)[^(k 2 )(jp 3 + m) ^(k 2 )(^' + m) ^{k^fa + m) / Al (fci)]u(p, s 



(p -J?i-fci)(p -P2 )(P -P3~ k 2 

In limit I, Eqs. ( ICip and ( lC4p can be added such that the denominator reduces to (p- ki) 2 (p 
k 2 ). Using Eqs. (lAIj) . (1A9| . flAlOl) and (HE}, we obtain 

[(5m 2 ) 3a + (5m 2 ) 3fe ] J = - / d 2 k 1± d 2 k 2± /" ^* 2+ 



2(2^) 6 y y 32^+^+ 

rr[ / ^ Al (A; 1 )(i) + m) / A2 (fc 2 )(^ 3 + m) jt X2 {k 2 ){£> + m) ^(fci)(^ + m)] 

(p • &i) 2 (p • fc 2 ) 

(C7) 
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After calculating the trace, Eq. ( IC7j) leads to 



[(Sm 2 ) 3a + (5m 2 ) 3b y 

-~Wfl dk ^ dk ^J If 14 m 

Again, in limit II, the denominator in the sum of Eqs. fICip and flU4p reduces to {p-ki)(p-k 2 ) 2 
leading to 

,2\ 1// e / j2i, j2i, / uhj l uh, 2 



[(5m% a + (dm^Y 1 = - —— / d 2 k 1± d 2 k 2± 



2(2^) 6 J J 32/fc+/c+ 

Tr[/ A2 (A; 2 )(^ + m) / Al A + m) 4 X2 (k 2 )(fa + m) / Al (fci)( j6 + m)] 



(p • fci)(p ■ 



2 



\2 



(C9) 



Calculating the trace, Eq. (1C9|) reduces to 

[(<W) 3a + (5m 2 ) 3b ] 7/ 

e4 ^ /" dkt [2(p ■ 6(A; 1 )) 2 (p ■ e(fc 2 )) 2 - (p ■ fcQ(p ■ e(k 2 )) 2 } 
= -j2^J dk ^ dk ^Jlflf 4HiPJ (C10) 

In limit III, we obtain 

e 4 /" ... ,0. f dkf dk+ (p-e(h)) 2 (p-e(k 2 )) 2 



!(*»*)». + (W = - w y <^u^ y ^"-^'^y (cm 

Similarly, T 3c leads to 
Sml =-J— R [ d 2 k 



^ f ^2i j2i /" dk^ dk 2 
3c =2(2^ 7 d kl±d k2± J 32kJkJpM? 



Tr[^(k 2 ){fa + m) jt x *(k 2 ){fa + m) 4 Xl {h)( fa + m) jt Xl (&i)(i> + m)] 



(C12) 



(p -Pi - ^ )(p - p' 2 )(p - p 3 - k 2 ) 
where pi and p 3 are defined by Eqs. ( 1C2I) and ( 1C5|) and p 2 = p. Note that this diagram 
is one-particle reducible, and therefore the energy denominator associated with the single- 



particle state vanishes. We shall use the Heitler method 
Using this method, we write 
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29[ for evaluating all such integrals. 



D 



29] 
1 



(p~ ~ P2 )(P~ ~Pi ~ k i )(P" -Ps ~ k 



2 



dp'-5(p'--p-)- ^— V — -, (C13) 

W - P2 ){P - Pi ~ k i ){P ~Ps -h) 



Using the relation between distributions 



? —Ap'~-P~) = -b'(p'--p~). (C14) 



{p'--p' 2 ) ^ 2 
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and integrating by parts we obtain 



dJ- J d j/- 5 tf--p-)JL 



1 



(p'~ - Pi - h ){v'~ -p 3 -k 2 

1 



2{p -Pi-^yip -p 3 -k 2 ) 2(p -p 1 -k 1 )(p -p 3 -k 2 ) 2 
Using Eq. O in Eq. flC15|> we obtain 



(C15) 



D 



Pi 



pi 



2(p-k 1 y( P -k 2 ) 2(p • h)( P ■ k 2 y 

Thus, in limit I, Eq. (1C12j) becomes 

dkf dk+ 2(p • e(A; 1 )) 2 (p • e(k 2 )) 2 - (p ■ k 2 )(p ■ e(h)) 2 



(C16) 



(27T) 



P 



+ 



8^ 



(C17) 



(p- h) 2 (p- k 2 ) (p ■ h)(p ■ k 2 ) 2 
Similarly, in limit II, Eq. ( 1C12I) leads to 

dkf dkt 2(p ■ e(k x ))\p ■ e(k 2 )f - (p ■ h){p ■ e(k 2 )) 2 



I3c\ 



(27T) 



^k 1± d Z k 2± 



Pi 



ip^ 



_(p ■ h) 2 (p ■ k 2 ) (p ■ h)(p ■ k 2 ) 2 
and in limit III, it gives 



(C18) 



i(Sm 2 h c } IH 



1 d 2 k 1± d 2 k 2± 1 dk ^ dk " 2 



(2tt) 6 J " " 1J -™ 7 fc+ kt 

4(p-fci)(p- k 2 ) 2 



4(p-A; 1 ) 2 (p- fca) 

The traces are calculated using Mathematica. 

The contribution corresponding diagrams in Fig. 4 is given by 



(C19) 



(<5m 2 ) 4 = (5m 2 ) ia + (<5m 2 ) 46 + (5m 2 



I ic 



In limit I, (<5m 2 ) 4a reduces to 
i(5m 2 ) 4a Y = 



r d 2 k l± d 2 k 2± f dk ? dk2 



(2tt)6 J "'^ ^ J k+ ^ 
[2p+(p ■ ejkMP ■ cihyMh) ■ e{k 2 )) + fcj> ■ 6(fc 1 )) 2 ] 
8p+(p • • k 2 ) 



(C20) 
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Note that (dm 2 )^ is not IR divergent when p ■ k\ 7^ even if p ■ k 2 = 0. 
In limit II, (Sm 2 )^ reduces to 



i(5m 2 ) 4b } n = JLj j d 2 k 1± d 2 k 2± j 



1 2 

[2pj(p ■ efaMp ■ e(k 2 ))(<h) • e(k 2 )) - kf(p ■ e(k 2 )) 2 ] 
8pf(p -ki)(p- k 2 ) 

Note that (dm 2 )^ is not IR divergent when p ■ k 2 7^ even if p ■ k± — 0. 

In limit III and after adding the contributions from Figs. 4(a) and 4(b), we get 



(C21) 



[(Sm% a + (5m% b } in = -JjL- j d 2 k 1± d 2 k 2± j 



dk+ dkj [(p ■ e(fci))(p ■ e(fc 2 ))(e(fci) ■ e(fc 2 ))] 
kf k+ Aip-kJip-h) 

(C22) 



For (Sm 2 )^, we use the Heitler method illustrated in Eqs. ( 1C13j) - flC15j) and obtain 



(Sm 2 ) 4c [ d 2 k 1± d 2 k 21 _ f dktdkt 



2(2tt) 6 J J 32k+k£p+ 

Tr[/P*{h x )(£ + m) ^(fci)(^ + m) ^ 2 (h)j+ / A2 (^)(i> + m)] 



(p • fo! 



2 



(C23) 



which finally leads to 



e 4 f , , , , f dki dk,2 {p ■ e{ki)) 2 



Similarly, contribution from Figs. 5(a)-(c) is given by 

(5m 2 ) 5 = (5m 2 ) 5(l + (<5m 2 ) 5b + (5m 2 ) 5c 
In limit I, (Sm 2 ) 5a reduces to 



l(5m% a Y = -J^e J d 2 k 1± d 2 k 2± J 



[2p+(p ■ 6(fcQ)(p ■ ejk^Mh) ■ e(fc 2 )) + fcj> ■ e(fcQ) 2 ] 
8p+(j9- fc 2 ) 

Similarly, in limit II we get 

[( W = - ^ / ,W 2 k 2 , / ff 

[2p+(p ■ 6(fc0)(p ■ e(fc 2 ))(e(A;i) ■ e(fr)) + fc+(p ■ e(k 2 )) 2 } 
8pf(p ■ kx)(p ■ k 2 ) 



(C25) 



(C26) 
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and in limit III we get 

e4 fj^ f dk t dkt [(p ■ e(A*))(p • e(fc 2 ))(e(A; 1 ) • e(k 2 ))] 



[(5m 2 ) 5a ] m = j^y 6 J d 2 k l± d 2 k 2± J 



kf fc+ 4( P .fc 1 )( P .fc 2 ) 



Taking limit I (5m 2 ) 56 reduces to 

p 4 f f rlh + rlh + 

,2\ 1/ e / j2i, r2i, / u,ft l u/i 2 



[(5m 2 ) 5 ^ = -J^y 6 J d 2 k 1± d 2 k 2± J 



"a ^2 

[2p+(p ■ 6(fc 1 ))(p ■ e(fc 2 ))(e(fcQ • e(fc 2 )) - k+jp ■ e(fc 1 )) 2 ] 
8^3 (P- h){p-k 2 ) 

Taking limit II we get 

&A f ,721, j2,, /" dkf dk 2 



i(5m 2 ) 5b ] n = ~J^y 6 J d 2 k 1± d 2 k 2± J 



[2pt(p ■ e(fci))(p ■ e(k 2 )Mh) ■ e{k 2 )) - fc+(p ■ £ (fc 2 )) 2 ] 
8pf(p -ki)(p- k 2 ) 

and taking limit III we get 

e 4 f J2l /■ dtf ^2 + [(P " e(*i))(p " • e(A*))] 



[(5m 2 ) 5fe ] m = - JL- I d 2 k 1± d 2 k 2± I 
(5m 2 ) 5 C in the limit I, reduces to 

Contribution of Figs. 6(a)-(c) is given by 

(5m 2 ) 6 = (5m 2 ) 6a + (5m 2 ) 6b + (<5m 2 ) 6c 
In limit I, (Sm 2 ) 6a reduces to 



K^U 1 = -j^y 6 J d 2 k l± d 2 k 2± J 



dk-^ dk^ 



1 2 

[2p + (p ■ e(h))(P ■ <k 2 )Mh) ■ e(k 2 )) + fffr ■ e(fcQ) 2 ] 
8p+(p • fci)(p • k 2 ) 

Note that (5m 2 )g a is not IR divergent when p ■ k\ 7^ even if p ■ k 2 = 0. 
In limit II, (5m 2 ) 6 f, reduces to 

,2\ 1// e4 /" j2i, >2i, /" dkf dk 2 



i(5m 2 U] H = - J d 2 k 1± d 2 k 2± J 



rti rug 

iM(P ■ <h)){p ■ e(k 2 )Mh) ■ e(fc 2 )) - kf(p ■ e(fc 2 )) 2 ] 
8p+(p ■ h)(p- k 2 ) 
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Note that (5m 2 )Qb is not IR divergent when p ■ k 2 7^ even if p ■ k\ = 0. 
In limit III, sum of ((5m 2 ) 6a and (Sm 2 )^ reduces to 



2\ Till 



[(5m ) 6a + (5m ) 66 ] 



(2tt)< 



(i 2 ki^(i 2 k2_L 



dk+ i(p ■ e(k l ))(p ■ e(k 2 ))(e(k 1 ) ■ e(k 2 ))} 



"'1 "'2 



4(p ■ ki)(p ■ k 2 ) 



For Fig. 6(c) we use Heitler method 28] to obtain the following result in limit I 



[(Sm 2 ) 6c Y 



(2k) 



d 2 \ii^d 2 ]i2± 



dkf dk+ (p-e(k 1 )f 
k+ k+ 8(p-k 1 ) 2 



(C34) 



(C35) 



Appendix D: Transition matrix in coherent state basis 

We will now present the calculation of (Sm 2 )'^ where prime denotes the extra contribu- 
tion arising due to use of coherent state basis. Contribution corresponding to Fig. 8 can be 
written as 

(5m 2 )' 8 = (5m 2 )' 8a + (Sm 2 )' 8b + (5m 2 )' 8c + (5m 2 )' 8d + (5m 2 )' 8e + (5m 2 )' 8f 

where 



16kfk£ 

Tr\^(k 2 )(fa + m) jt^(k 1 )(i) 1 + m) / Al (&i)(i> + m)](p ■ e Al (h))Q A (k 



(p-k 1 ) 2 [(p-k 1 ) + (p-k 2 )-(k 1 -k 2 )] 



(Dl) 



v 786 2(2^) 6 J J lQk+k+ 



Tr[^(k 2 )(fa + m) 4 Xl (ki)(h + m) ^(k 2 )(j> + m)](p ■ e x ^(k l ))Q A (k l ) 

(p-k l )(p-k 2 )[(p-k 1 ) + (p-k 2 )-(k l -k 2 )} 1 ] 

Calculating the trace, adding Eqs. (ID 1 [) and (1D2I) and taking limit I, we obtain Eqs. fl85|) . 
Similarly on can obtain the expressions for Figs. 8(c)-8(f) in appropriate limits. Fig. 9 
corresponds to the transition matrix element T 9 and its contribution is given by 

(5m 2 )' 9 = (5m 2 )' 9a + (6m 2 )' 9b + (5m 2 )' 9c + (5m 2 )' 9d 

where 

r^y im_ e4 /" ,2. ,2, / dkf dkj (p ■ e(k l )) 2 (p ■ e(k 2 )) 2 
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d 2 k 1± d 2 k 2± 



dkjdkt (p-e(fci)) 2 (p-e(fc 2 )) ; 
kf k+ Mp-k x )\p-k 2 ) 



e A (fci)e A (A: 2 ) (D4) 



2\' 



[(5m% c + ( Sm% d ] 



d 2 ki^d 2 k 21 ^ 



(p-e(h)) 2 (p-e(h)) 



dk^ dk 2 



(p ■ ejk^fjp ■ ejk^Y 
4(p- h) 2 (p- k 2 ) 



+ 



©A(fci)e A (fc 2 ) 



(D5) 



4(p- /c 2 ) 2 
Here we have used Heitler method to get the above result. 

Similarly, contribution coming from diagrams in Figs. 10 and 11 can be easily evaluated 

by putting the appropriate limits. 
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